function fent = drive(inparm)
 global count fout
n = length(inparm);
inparm = inparm';
count = count + 1;
%fent = -nlev_unew(n,inparm)
fent = -twoqu_u(n,inparm)
%fent = -nlev_u(n,inparm);
%fent = nlev_0(n,inparm)
save('outfile','fent','count','inparm')
fout(count) = fent;
%save('infile','n','inparm','-ascii')
%nlev_0;
%system('a.out');
%out = -load('fort.666')
%nput('');

%if mod(count,10)==0
 %   load('fort.111')
 %   figure(1)
 %   plot(fort(:,1),fort(:,5))
 %   hold on
 %   figure(2)
 %   plot(fort(:,1),fort(:,4))
 %   hold on
 %   figure(3)
 %   plot(fort(:,1),fort(:,3))
 %   hold on
%end
end

